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1.  Introduction 

We  axe  concerned  with  hyperbolic  systems  of  conservation  laws  with  umbilic 
points.  A  point  in  state  space  Uq  €  Rn  is  called  an  umbilic  point  for  a  system 
of  hyperbolic  conservation  laws 

(1.1)  dtU  +  dxF{U)  =  0,  Ue  R", 

if  some  wave  speeds  coincide  at  this  point  I7o,  that  is,  at  least  two  eigenvalues 
\i(U)  and  \j(U),i  ^  j,  among  the  n  real  eigenvalues  of  the  matrix  VF(? 7)  axe 
such  that  \i(Uo)  =  Xj(Uo)-  Such  umbilic  points  allow  a  degree  of  interaction, 
or  nonlinear  resonance,  between  distinct  modes,  and  lead  to  high  singularities, 
which  is  missing  in  the  strictly  hyperbolic  case. 

The  study  of  singular  hyperbolic  equations  has  an  extensive  history  dating 
back  at  least  to  the  work  of  Euler  [Eu]  two  hundred  years  ago,  when  he  proposed 
a  singular  equation  —  the  Euler-Poisson-Darboux  equation.  Its  close  relations 
with  wave  theory,  fluid  dynamics  as  well  as  geometry  have  attracted  great  atten¬ 
tions  from  mathematicians  for  two  centuries  including  Poisson  (1823),  Darboux 
(1914),  Riemann  (1860),  Volterra  (1892),  Beltrami  (1880)  and  also  Weinstein, 
Erdelyi,  Lions  and  others  in  »,he  60’s  (see  [Y]).  Recently,  a  study  of  the  behavior 
of  its  solutions  led  to  a  solution  of  the  nonlinear  system  of  isentropic  gas  dynam¬ 
ics  [Chi,  DCL]  (also  see  [Di2,  Ch2])  because  of  its  close  relation  with  entropy 
information.  On  the  other  hand,  the  theory  of  linear  equations  with  multiple 
characteristics  is  also  well  developed.  An  important  feature  of  such  equations  is 
the  loss  of  differentiability  [DeG],  which  leads  to  ill-posedness  in  Sobolev  spaces 
and  to  the  use  of  Gevrey  Classes  [Ge].  Another  feature  is  that  sign  conditions 
on  the  subprincipal  symbol  play  an  important  role  (cf.  [Fr,  Oh,  H]). 

Recently,  nonlinear  hyperbolic  systems  with  umbilic  points  have  arisen  from 
such  areas  as  multiphase  flows  in  porous  media,  elasticity,  water  wave  problems, 
and  maganetohydrodynamics.  Such  umbilic  points  appear  naturally  in  multi¬ 
dimensional  systems  of  conservation  laws.  In  particular,  Lax  showed  that  in 
three  space  variables  there  must  be  umbilic  points  if  the  number  of  equations 
is  2 (mod  4)  [Lai].  The  theory  of  local  solutions  for  such  systems  is  well  devel¬ 
oped  because  many  tools  for  linear  equations  still  can  be  used.  However,  since 
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such  systems  axe  nonlinear,  even  if  the  initial  data  are  smooth,  the  solutions  of 
the  Cauchy  problem  generally  develop  singularities  and  become  discontinuous  in 
finite  time.  This  is  an  exact  reflection  of  the  physical  phenomena  of  the  break¬ 
ing  of  waves  and  the  development  of  shock  waves.  An  effort  has  been  made  to 
understand  the  Riemann  solutions  for  such  systems  (cf.  [Gl,  IMPT,  IT,  SS2, 
SSMP]).  Two  kinds  of  degeneracy  are  classified,  which  govern  different  behavior 
of  solutions  near  umbilic  points.  A  typical  example  of  parabolic  degeneracy  is 
the  system  of  isentropic  gas  dynamics  (cf.  [Ch2]).  The  most  simple  example 
of  hyperbolic  degeneracy  is  the  system  with  a  rotational  symmetry,  one  linear 
degenerate  characteristic  field,  and  one  contact  characteristic  field  (cf.  [KK,  LW, 
F,  Ch3]). 

This  is  the  first  in  a  series  of  our  papers.  In  this  series  of  papers  we  focus 
on  isolated  umbilic  points  with  hyperbolic  degeneracy.  Near  such  an  isolated 
umbilic  point  one  can  scale  and  blow  up  singularity  to  yield  generically  a  homo¬ 
geneous  polynomial  flux,  determined  by  the  lowest-order  nontrivial  terms,  by  a 
Galilean  transformation.  For  the  2x2  case,  this  process  generically  leads  to  a 
homogeneous  quadratic  polynomial  flux.  Such  a  polynomial  flux  contains  some 
inessential  scaling  parameters  and  the  selection  of  a  unique  flux  from  each  equiv¬ 
alence  class  is  the  problem  of  normal  forms.  It  was  solved  by  Issacson,  Plohr, 
and  Temple  and  in  a  more  satisfactory  form  by  Schaeffer  and  Shearer  [SSI].  The 
classification  of  the  geometry  of  rarefaction  curves  and  some  preliminary  tools 
for  the  analysis  of  Riemann  problems  for  the  quadratic  flux  systems  are  also  pre¬ 
sented  in  [SSI].  The  Riemann  solutions  for  such  systems  were  constructed  by 
Issacson,  Marchesin,  Paes-Leme,  Plohr,  Schaeffer,  Shearer,  Temple  and  others 
(cf.  [IMPT,  IT,  SS2,  SSMP]). 

The  global  existence  of  weak  solutions  to  the  Cauchy  problem  for  a  special 
case  of  such  quadratic  flux  systems  was  solved  in  [K]  using  the  viscosity  method. 
The  method  of  compensated  compactness  was  used  and  the  main  analysis  in¬ 
volves  detailed  estimates  and  characterization  of  the  singularities  of  solutions  to 
the  associated  entropy  equation  which  is  of  Euler- Poisson-Darboux  type.  Classes 
of  Goursat  data  were  then  carefully  chosen  to  cancel  these  singularities  in  the 
construction  of  general  classes  of  regular  entropies.  In  [Rb],  the  analysis  in  [K] 
was  applied  to  a  slightly  different  system.  And  in  [Lu],  a  different  proof  was 
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given  independently  to  the  same  problem  studied  in  [K]. 

In  this  first  paper  of  this  series,  we  establish  an  L°°  compactness  framework 
for  sequences  of  approximate  solutions  to  general  nonstrictly  hyperbolic  systems 
with  umbilic  points.  Some  techniques  in  compensated  compactness  [Ta,  Mu,  Se] 
are  used  and  the  analysis  of  singularities  in  [K]  is  generalized  to  achieve  this 
framework.  Under  this  framework  approximate  solution  sequences,  which  are 
apriori  bounded  in  L°°  and  which  produce  correct  entropy  dissipations,  lead  to 
the  compactness  of  the  corresponding  Riemann  invariant  sequences.  In  another 
series  of  papers,  we  are  going  to  develop  the  corresponding  Lp  theory.  One 
of  the  principal  difficulties  associated  with  such  systems  is  the  general  lack  of 
enough  classes  of  entropy  functions  that  can  be  verified  to  satisfy  certain  weak 
compactness  conditions  in  the  div-curl  lemma  of  Tartar  [Ta]  and  Murat  [Mu]. 
This  is  due  to  possible  singularities  of  entropy  functions  near  the  regions  of 
nonstrictly  hyperbolicity.  The  analysis  leading  to  the  compactness  involves  two 
major  steps: 

In  the  first  step,  we  construct  regular  entropy  functions  governed  by  a  highly 
singular  entropy  equation.  There  are  two  main  difficulties.  The  first  is  that,  in 
general,  the  coefficients  of  the  entropy  equation  are  multiple-valued  functions 
near  the  umbilic  points  in  the  Riemann  invariant  coordinates,  which  is  missing 
in  the  special  cases  [K,  Rb].  This  difficulty  is  overcome  by  a  detailed  analysis  of 
the  singularities  of  the  Riemann  function  of  the  entropy  equation  in  Section  3 
and  Section  4.  This  analysis  involves  a  study  of  a  corresponding  Euler- Poisson- 
Darboux  equation  using  a  majorant  idea  and  requires  very  complicated  estimates 
and  calculations.  Finally,  an  appropriate  choice  of  Goursat  data  leads  to  a 
cancellation  of  singularities  and  we  obtain  regular  entropies  in  the  Riemann 
invariant  coordinates.  The  second  difficulty  is  that  the  nonlinear  correspondence 
between  the  U -coordinates  and  the  Riemann  invariant  coordinates  is,  in  general, 
irregular.  A  regular  entropy  function  in  the  Riemann  invariant  coordinates  is 
usually  no  longer  regular  in  the  physical  coordinates  U.  We  overcome  this  by  a 
detailed  analysis  of  the  correspondence  between  these  two  coordinates. 

In  the  second  step,  we  study  the  structure  of  the  Young  measures  associ¬ 
ated  with  the  approximate  sequences,  and  prove  that  the  support  of  the  Young 
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measures  lies  in  finite  isolated  points  or  separate  lines  in  the  Riemann  invariant 
space.  This  is  achieved  by  a  delicate  use  of  Serre’s  technique  [Se]  and  regular 
entropy  functions,  constructed  in  the  first  step,  in  the  Tartar-Murat  commuta¬ 
tion  equation  [Ta]  for  Young  measures  associated  with  the  approximate  solution 
sequences. 

This  compactness  framework  is  successfully  applied  to  prove  the  conver¬ 
gence  of  the  Lax- Friedrichs  scheme  [La3],  the  Godunov  scheme  [Go]  and  the 
viscosity  method  for  a  canonical  class  of  the  quadratic  flux  systems  and  other 
related  systems  in  Section  7.  Corresponding  existence  theorems  of  global  en¬ 
tropy  solutions  for  such  systems  are  established.  The  compactness  is  achieved 
by  reducing  the  support  of  the  corresponding  Young  measures  to  a  Dirac  mass 
in  the  physical  space.  This  involves  a  careful  construction  of  many  other  special 
solutions  to  the  nonstrictly  hyperbolic  entropy  equations  and  a  detailed  analysis 
of  finite  Tartar-Murat  functional  equations. 

In  the  forthcoming  second  paper  [CK1]  we  will  apply  this  framework  to 
solve  the  other  three  remaining  canonical  classes  of  quadratic  flux  systems  and 
related  systems.  Several  convergence  theorems  of  L°°  approximate  solutions  and 
existence  theorems  of  global  solutions  will  be  presented. 

In  the  third  paper  of  this  series  [CK2],  we  will  solve  the  system  of  three- 
phase  incompressible  flows  in  porous  media,  studied  by  many  authors,  with  the 
aid  of  a  variant  of  this  framework.  By  a  detailed  analysis  of  the  geometry  of 
wave  curves  and  the  dissipation  of  entropy  waves  we  obtain  the  convergence  of 
approximate  solutions  (such  as  the  Lax- Friedrichs  scheme,  the  Godunov  scheme, 
and  the  viscosity  method)  and  the  existence  of  globed  entropy  solutions  for  the 
Cauchy  problem  with  large  initial  data  for  this  physical  system. 

In  connection  with  earlier  work  on  compactness  frameworks  on  approximate 
solutions  to  strictly  hyperbolic  conservation  laws  we  refer  the  reader  to  the  work 
of  Tartar  [Ta]  for  scalar  conservation  laws,  and  to  the  work  of  DiPerna  [Dil]  and 
Serre  [Se]  for  2  x  2  systems.  DiPema’s  analysis  [Dil]  is  based  on  a  study  of 
the  Lax  progressing  entropy  waves  in  state  space  [Lai],  and  in  particular,  on 
relationships  between  their  structure  and  the  structure  of  the  Young  measures 
by  using  the  strict  hyperbolicity  and  the  convexity  of  the  systems.  Serre  [Se] 


provided  another  approach  to  establish  such  a  compactness  framework  for  gen¬ 
eral  systems  by  using  certain  kinds  of  Goursat  entropies  whose  regularities  are 
ensured  by  the  strict  hyperbolicity  of  the  systems. 

Regarding  work  on  other  related  topics  on  hyperbolic  conservation  laws  with 
umbilic  points  we  refer  the  reader  to  Isaacson,  Maxchesin,  and  Plohr  [IMP]  for 
a  general  description  on  transitional  shock  waves,  to  Liu  [Liu]  for  a  discussion 
on  asymptotic  stability  of  such  waves,  and  to  Glimm  [Gl,  G3]  for  a  review 
of  applications  of  bifurcation  theory  and  geometry  to  the  analysis  of  Riemann 
solutions. 


2.  The  Classification  of  Quadratic  Fluxes  and  Wave  Curves 

Consider  a  hyperbolic  system  of  conservation  laws 

(2.1)  dtU  +  dxF(U)  =  0,  U  €  R2, 

with  an  isolated  umbilic  point.  An  point  Uo  €  R2  is  called  an  isolated  umbilic 
point  if  VF(Uo)  is  diagonalizable,  and  there  is  a  neighborhood  N  of  Uq  such 
that  VFt(U)  has  distinct  eigenvalues  for  all  U  €  N  —  Uo,  where 

VFt(U)  =  F(Uo)  +  VF(Uo)(U  -  Uo)  +  l-(U  -  U0)TV2F(U0)(U  -  U0). 

Take  the  Taylor  expansion  for  F{U)  about  U  =  Uq\ 

(2.2)  F(U)  =  F(U0)  +  VF(U0)(U  -UQ)  +  hu-  U0)T' V2F(U0)(U  -U0)  +  h.o.t. 

where  h.o.t.  represents  the  remainder.  The  flux  function 

(2.3)  Q(U  -  Uo)  =  F{Uq)  +  VF(U0)(U  -  U0)  +  \{U  -  U0)r V2 F(U0)(U  -  U0) 

determines  the  local  behavior  of  hyperbolic  singularity  near  the  umbilic  point 
Uo.  Since  VF(Uo)  is  diagonalizable,  we  can  make  a  coordinate  transformation 
to  eliminate  the  linear  term  from  (2.3)  and  relabel  U  —  Uo  as  U  to  obtain 


(2.4) 


dtU  +  dxQ{U)  =  0, 
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from  (2.1)  and  (2.3),  where  Q{U )  =  \UT V2 F(Uo)U ’.  From  the  normal  form 
theorem  in  [SSI],  there  is  a  nonsingular  linear  coordinate  transformation  to 
transform  the  system  (2.4)  into 

(2.5)  dtU  +  dx(dC(U))  =  0, 
where 

C(U)  =  -(-auz  +  bu2v  +  uv2),  a  ^  1  +  b2 . 

A  o 

We  now  analyze  the  geometry  of  rarefaction  wave  curves  and  the  genuine 
nonlinearity  of  the  quadratic  flux  system  (2.5)  in  the  sense  of  Lax  [La4].  We  will 
make  use  of  this  information  to  understand  the  classification  of  the  canonical 
form  into  four  regions  in  the  (a,  6)-plane  and  thereby  rederiving  some  of  the 
results  in  [SSI]  from  a  different  viewpoint. 

This  analysis  is  also  essential  in  obtaining  an  L°°  apriori  estimate  for  the 
sequences  of  viscous  and  finite  difference  approximate  solutions  via  invariant 
region  techniques.  We  establish  such  estimates  in  Section  7.  For  the  symmetric 
case  of  (2.5)  ( b  =  0),  the  existence  of  invariant  regions  and  the  geometry  of  wave 
curves  were  analyzed  in  [K].  We  remark  here  that  our  analysis  will  be  consistent 
with  Darboux’s  local  analysis  near  an  isolated  umbilic  point  (see  [Da]). 

Recall  from  (2.5)  that  the  flux  vector  and  matrix  take  the  form: 

(2.6)  F(U)  =  dC(U)  =  ^(au2  +  2  buv  +  v2,bu 2  +  2uu)T, 
and 

(2.7)  W)-(“t  1  ‘T 

The  eigenvalues  and  eigenvectors  Eire  from  (2.7), 

(2.8)  A,  =  ^[(a  +  l)u  +  bv  +  (  —  1)* \/ ((a  -  l)it  +  bv )2  +  4(6it  +  u)2],  i  =  1,2. 
and 


(2.9)  r,  =  ((a  -  l)u  +  bv  +  (-1)'  \/((a  -  l)u  +  bv)2  +  4(6u  +  v )2  ,  2(bu  +  v))T , 
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respectively.  It  is  then  immediate  that  as  long  as  a  ^  1  +  b2,  Aj  =  A2  *£=> 
( u,v )  =  (0,0),  so  that  (0,0)  is  the  unique  umbilic  point  for  (2.5). 

The  jth  family  of  rarefaction  curves  Rj  is  defined  as  the  family  of  integral 
curves  of  the  vector  field  given  by  rr  Therefore,  Rj  is  defined  by  the  following 
ordinary  differential  equation: 


du  (a  —  1  )u  +  bv  +  (— 1  y  y/((a  —  l)u  +  bv )2  +  4(6u  +  v )2 
dv  2  (bu  +  u) 

—Fj(u,v) 

(2-10) 

(a  —  l)o:  -f  b  +  (  — 1  y  sign(v)y/[(a  —  1)  +  ba]2  +  4(6  +  a)2 


where  a  =  ^.  Define 


(2.11) 


-  l)a  +  b]  +  (  — 1)J  \/[{a  -  1)  +  ba]2  +  4 (b  +  a)2 
2(ba  +  1) 


Then  we  have 


(2.12) 


v  >  0, 
v  <  0. 


To  analyze  the  geometry  of  the  Rj  curves,  we  start  by  noting  that 


(2.13) 

du 

=  Fj(u,v), 

uV 

<Pu 

d2v 

=  -daFj(u 

V 

(2.14) 

=  -daFj(u 

V 

,  du 
'  dv 


when  v  ^  0. 


We  collect,  in  the  following  lemma,  some  basic  properties  of  g}. 
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Lemma  2.1.  The  function  gj(a),j  =  1,2,  satisfies 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


Next,  we  make  use  of  (2.11)  —  (2.14)  and  Lemma  2.1  to  obtain  a  qualitative 
picture  of  the  R,  curves.  We  distinguish  several  cases.  We  remark  that  from 
(4)-(6)  in  lemma  2.1,  a  =  1  +  62  defines  a  boundary  curve  in  the  (a,  6)-plane 
separating  qualitatively  different  wave  curve  geometries.  We  shall  elaborate  on 
this  below. 

First,  we  assume  that  a  >  1  +  62.  Using  Lemma  2.1,  we  obtain  the  following 
graphs  of  g}  in  the  two  subcases  6  >  0,  and  b  <  0. 

(note  to  printer:  place  Fig.  2.1  here  in  text) 

For  simplicity,  we  first  consider  the  case  b  >  0.  From  (2.14),  the  change  in 
convexity  (as  a  function  of  the  slope)  of  the  wave  curves  depends  on  the  location 
and  number  of  roots  of  gj (a)  —  a.  A  computation  using  the  formulae  for  the 
g'jS  shows  that,  in  the  present  case,  these  locations  are  given  by  the  roots  of  the 
cubic  polynomial 


h(a)  =  —6a3  4-  (a  —  2)a2  +  26a  +  1. 
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The  discriminant  of  h  is  given  by  A  =  —  32£>4  +  b2( 27  +  36(a  —  2)  —  4 (a  — 
2)2)  +  4(a  —  2)3.  Thus,  A  =  0  gives  a  new  boundary  in  the  (a,  6)-plane  which 
distinguishes  different  wave  curve  geometries.  This  corresponds  to  the  division 
between  regions  III  and  IV  in  [SSI].  When  A  <  0,  h  has  three  real  roots  ao,«i, 
and  c*2 .  We  therefore  obtain  the  following  figures  for  the  Ry  curves  in  the  case 
b  >  0,a  >  1  +  b2,  A  <  0. 

(note  to  printer:  place  Fig.  2.2  here  in  text) 

The  case  b  <  0,a  >  1  -f-  i2,  A  <  Ois  completely  similar.  We  show  the 
corresponding  graphs  below  in  Fig.  2.3. 

(note  to  printer:  place  Fig.  2.3  here  in  text) 

When  A  >  0,  h  has  only  one  real  root  and  the  wedge-shaped  regions  in  Fig. 
2.2  and  Fig.  2.3  collapse  into  a  line.  This  corresponds  to  region  IV  in  [SSl]. 

(note  to  printer:  place  Fig.  2.4  here  in  text) 

Next,  we  consider  the  case  a  <  1  +  b2.  Using  Lemma  (2.1),  the  gj  diagrams 
are  as  follows: 

(note  to  printer:  place  Fig.  2.5  here  in  text) 

By  completely  imilar  reasoning  as  in  previous  cases,  we  obtain  the  following 
R}  curves  diagrams: 

(note  to  printer:  place  Fig.  2.6  here  in  text) 

We  now  turn  to  investigate  the  genuine  nonlinearity  in  the  sense  of  Lax 
[La4]  for  the  quadratic  system  (2.5).  It  will  turn  out  that  genuine  nonlinearity 
allows  a  breakdown  of  the  last  case  above  (a  <  1  -f  b2)  into  two  subcases.  By 
definition,  (2.5)  is  genuinely  nonlinear  in  the  jth  characteristic  field  at  a  point 
(u,  v )  if  Tj  •  VA j  0,  i  t ^  j,  at  ( u ,  u).  A  calculation  using  (2.8)  —  (2.9)  shows  that 

rrV\j  =aC  +  Zby  +  (-iy^^ 

C  =  (a  —  l)u  +  bv, 

(2.16)  y  =  bu  +  v. 


36yC+2(a  +  3)t/2 
V/C2+4y2 
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Using  (2.16),  we  find  that 


(2.17) 


r,  •  V\, 


y  =  o,(— i)-*C  <  o, 


'  y  =  0,j,=  5i±=2c, 

,y  =  0,y=  'K 


a  >  I b2, 
a=|  &2, 
C,  a  <  |62. 


where  D  =  —  12(a  —  |62)(a  -j-  3)2.  (2.17)  shows  that  the  curve  a  =  |62  divides 
the  region  {(a,  6)|a  <  1  +  b2}  into  two  subregions  according  to  a  global  change 
in  loci  of  loss  of  genuine  nonlinearity.  This  corresponds  to  the  division  between 
region  I  and  II  in  [SSlj. 


We  summarize  the  four  boundaries  in  the  (a,6)-plane  separating  different 
behavior  of  the  wave  curves  and  the  qualitative  change  in  the  system  (2.5)  as 
one  crosses  these  boundaries  in  the  following  table: 

(note  to  printer:  insert  Table  2.1  here  in  text) 


(note  to  printer:  insert  Fig  2.7:  Four  Boundary  Curves  here  in  text) 


3.  Riemann  Invariants  and  Genuinely  Nonlinearity  for  the  Quadratic 
Flux  System 

In  this  section,  we  study  the  Riemann  invariants  of  the  quadratic  flux  system 
(2.5).  We  also  study  the  monotonicity  of  A j,i  =  1,2,  as  a  function  of  Riemann 
invariants  Wj,j  =  1,2.  We  remark  that  this  does  not  follow  from  our  knowledge 
of  genuine  nonlinearity  (Section  2)  as  the  map  J  :  (u,v)  — ►  ( wi,wi )  is  neither 
C1  nor  globally  invertible  in  general. 

Riemann  invariants  Wj  —  Wj(u,v),j  =  1,2,  are  defined  as  functions  that 
are  constants  along  any  rarefaction  wave  curves  of  the  ith  family  R,  where  z  /  j. 
On  regions  where  w}  is  differentiable,  it  is  easy  to  check  that  since  R;  curves 
are  integral  curves  of  the  vector  field  r^,  iq  •  Vu?;  =  0,  i  ^  j. 

We  shall  use  the  ordinary  differential  equation  (2.10)  to  define  w}.  Care 
must  be  taken  in  its  definition  to  make  sure  that  it  comes  out  as  a  single- valued 
function  globally.  For  simplicity,  we  restrict  ourselves  to  a  half  plane  domain 
defined  below.  This  domain  is  also  an  invariant  region  for  the  viscous  system 
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associated  with  (2.5)  (see  Subsection  7.1).  We  renark  that  the  w'jS  are  not 
uniquely  defined.  We  shall  make  the  choice  that  guarantees  maximal  regularity 
of  certain  quantities.  We  elaborate  on  this  below.  The  advantage  of  this  choice 
will  become  clear  in  Section  4  when  we  study  the  entropy  functions.  From  now 
on,  we  assume  that  6  >  0  for  simplicity.  Another  case  is  completely  similar. 

Proposition  3.1.  Consider  (2.5)  in  region  TV,  i.e.,  A  >  0  .In  this 
case,  the  cubic  polynomial  h(a )  has  only  one  real  root  a0.  Denote 

.  j3  4-  46o:o  +  (a  —  2)a;2| 

iao!\/((a  ~  l)ao  +  b)2  -f-  4(fto  +  l)2 

Consider  the  half  plane  domain  lk  —  {  ( u ,  u)  )  (  —  l)fe(u  —  oo^’)  >  0  },  k  = 

1, 2.  Then  the  following  formulae  define  a  pair  of  Riemann  invariants 
for  (2.5)  on  Tk: 

(3.1.1)  Wj(u,v)  =  f  Hi(a)  da},  i^j, 

1  +  J  0 

with 


(3.1.2) 


Ht{d)  = 


(1  +  a2)  E{&) 
2(ao  +  d)  D(a)  ’ 


where 


E(a)  =  —  2 b(a0d  —  l)2  +  (a  —  3)(ao  -f  d)(aod  —  1)  +  b(a0  +  d)2 
+  (-l)«+1(a0  +  d)v/Q(d), 

Q(a)  =[(a  -  l)(a0d  -  1)  +  b(a0  +  d)]2  +  4[6(a0d  -  1)  +  (a0  +  d)]2, 
D(d)  =  -  fe(a0d  -  l)3  +  (a  -  2)(ao  +  d)(a0d  -  l)2 

+  26(q0  +  d)2(a0d  -  1)  +  (a0  +  a)3, 

1 

u  =u  H - v, 

1 

V  =v - u , 

<*0 

u 

a  =—. 
v 
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Next,  we  study  the  ratio  ^  asa  function  of  a  =  jr.  We  show  that  our  choice 
of  the  definition  of  uj' s  are  natural  in  the  sense  that  such  Riemann  invariants  are 
globally  well-defined  on  2*  and  this  ratio  is  an  real  analytic  function  of  <5  in  a 
maximal  domain.  This  analyticity  is  crucial  in  analyzing  the  entropy  equations 
of  (2.5)  and  is  indispensable  in  the  proof  of  some  of  our  main  results  (theorem 
4.4.3  and  4.5.2). 


Proposition  3.2.  Let  iv\  and  w2  be  defined  as  in  Proposition  3.1 
Then 

—  =  r  j(a,  sgn(v)) 

Wi 


(3.1.3) 


H  ~exi){(—iy+1 0  r  i+j,  ij  =  1,2, 

Jo  L>(a) 

where  the  functions  Ti  and  r2  are  real  analytic  in  a  €  R’J  {(— l)*oo}. 


The  next  proposition  concerns  the  monotonicity  of  the  wave  speed  A,  in 
the  variable  w,.  This  is  important  in  the  reduction  of  the  Young  measures  for 
approximate  solution  sequences  in  Section  7. 


Proposition  3.3.  Suppose  that  A  >  0.  Given  v’i  and  tu2  as  defined  in 
Proposition  3.1,  the  eigenvalues  X,,i  =  1,2,  are  well-defined  functions 
in  (wi,iv 2)  €  J(Ik)-  Moreover,  we  have 
d\ 

0^-  #  0,  for  all  ( wx,w2 )  €  J(Tk)  -fai  =  ^  i  =  1,2. 


4.  Entropy  Functions 
4.1.  The  Entropy  Equation 

Following  the  standard  definition  in  [La2],  we  call  a  pair  of  scalar  functions 
( r](u,v ),  q(u,v ))  an  entropy-entropy  flux  pair  for  (2.1)  if  for  smooth  solutions 
( u(x ,  t),  v(x,  t ))  of  (2.1),  we  have  the  additional  conservation  law 

dtr](u(x,t),v(x,t))  +  dzq(u(x,t),v(x,t))  =  0. 
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It  is  easy  to  check  that  this  happens  iff  r?  and  q  satisfy  the  compatibility  condition 
(4.1.1)  V77VF  =  Vq. 


Eliminating  q ,  we  get  a  second  order  equation,  the  entropy  equation, 


(4.1.2)  Qu^Juu  4"  (<7t>  fu)T}uv  fv'Hvv  4“  (juti  fu t))^?a  4"  {guv  fvv^v  —  O5 
where  we  denote  F{u,v)  =  (f(u,v),g(u,v))T. 


(4.1.2)  is  a  linear  hyperbolic  equation  whose  characteristic  variables  turn 
out  to  be  the  Riemann  invariants.  A  simple  calculation  gives  the  characteristic 
form  as 


(4.1.3) 


f]wi  U>2 


4- 


A2  -  A,  ’lm 


4.2.  The  Entropy  Equation  for  the  Quadratic  Flux  System 

We  now  investigate  the  properties  of  the  coefficients  of  (4.1.3)  in  the  case 
of  the  quadratic  flux  system  (2.5).  It  will  turn  out  that  these  coefficients  cannot 
be  written  down  explicitly  as  functions  of  W\  and  w2  in  any  reasonable  closed 
form.  Instead,  we  will  compute  them  as  functions  of  a  =  j  and  study  their 
properties  using  information  from  Section  3.  We  summarize  the  main  results  in 
the  following  proposition: 

Proposition  4.2.1.  Consider  the  quadratic  dux  system  (2.5).  Sup¬ 
pose  that  we  axe  in  region  TV  in  the  (a,  b)-plane  so  that  A  >  0.  Then 
the  coefficients  of  (4.1.3)  satisfy,  for  i  ^  j, 

(1) 

*jwi  _  Tj(d,sign{v )) 

A2  —  Ai  w, 

Ti(a,  sign(v))  =  -Tj(d,  -sign(v)), 

where  Tj(q,±1)  and  Tj(a,±l)Tj(a,±l)  are  real  analytic  in  a  € 

R  U  {±00}; 
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(2)  If  v  >  0,  then 

Tj(a,i'zn(a)) 

Wi  W2  ~  U)\ 


W2  —  Wi 


Here,  Aj(cr)  =  Aj(^),  j  =  1,2,  are  real  analytic  in  cr  €  RU  {±00} 
and  i  €  RU  {±oo},  respectively; 

(3)  There  exists  M  >  0  such  that  ||A(^)||oo,  and 

||(wi,  w2)V2A(^-)(^i,  w2)t||oo,  i,j  =  1,2,  are  all  bounded  by  M 
for  all  w\  and  u;2. 

The  bounds  on  Ai  in  (3)  of  Proposition  4.2.1  will  prove  to  be  crucial  in 
analyzing  the  singular  behavior  of  77.  We  shall  elaborate  on  this  in  Section  4. 
As  a  corollary  of  (2)  above,  we  obtain 

Lemma  4.2.2.  The  entropy  equation  for  the  quadratic  flux  system 

(2.5)  takes  the  form 

.  M%)  M%) 

(4-1.4)  T]WlWi  +  - —Vwt - —Tlrv,  =  0. 

W2  —W 1  W2  ~  W\ 

It  is  clear  from  (4.1.4)  that  the  coefficients  of  the  entropy  equation  become 
singular  along  the  line  w  1  =  wj.  For  the  quadratic  flux  system,  we  have  the 
relation  uq  <  0  <  u;2  (see  Section  3).  Therefore,  upon  restriction  to  the  physical 
domain,  the  coefficients  are  singular  only  at  the  umbilic  point.  The  singular  na¬ 
ture  of  the  coefficients  also  suggest  that  general  solutions  77  cannot  be  expected 
to  be  C2.  This  presents  major  difficulties  in  verifying  H~l  compactness  condi¬ 
tions  in  the  div-curl  lemma.  This  is  reminiscent  of  similar  difficulties  for  proving 
convergence  of  approximate  solutions  to  the  gas  dynamics  equations  (see  [Ch2, 
Di2] ).  However,  the  singularities  of  the  entropies  are  of  somewhat  different  char¬ 
acters  in  the  two  cases  and  the  methods  for  resolving  these  difficulties  are  not 
the  same.  The  method  we  shall  present  generalizes  that  in  [K]  and  consists  of 
constructions  of  very  general  classes  of  nonsingular  (C2)  entropies.  This  anal¬ 
ysis  is  not  restricted  to  the  quadratic  flux  systems  and  can  be  generalized  to 
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obtain  a  compactness  framework  theorem  (Theorem  6.3)  for  convergence  of  ap¬ 
proximate  solutions  to  general  nonstrictly  hyperbolic  systems  with  an  isolated 
umbilic  point. 

We  remark  that  (2)  in  Proposition  4.2.1  implies  that  is  real  analytic 

in  u>i  and  iu2  except  at  the  umbilic  point  (tci,u>2)  =  (0,0)  and  that  Ai(^)  is 
multi-valued  at  the  umbilic  point.  We  now  interpret  this  by  comparing  (4.1.4) 
to  the  classical  Euler- Poisson-Darboux  equation.  The  classical  EPD  equation  in 
characteristic  form  is 

(4.1.5)  T)w\wi  "h  •Puii  =  0. 

W2  ~  W i  W2  ~  W\ 

where  j3\  and  /32  are  constants.  The  significance  of  these  constants  he  in  the 
fact  that  they  completely  determine  the  singular  behavior  of  solutions  to  (4.1.5). 
The  EPD  equation  arises  as  the  entropy  equation  for  the  isentropic  gas  dynamics 
equations  (see  [Di2,  Ch2]).  A  comparison  between  (4.1.4)  and  (4.1.5)  indicates 
that,  heuristically,  we  should  expect  the  singularity  of  the  entropy  near  the 
umbilic  point  depends  on  the  angle  of  approach  and  the  “size”  of  At .  This  turns 
out  to  be  the  correct  picture  and  motivates  much  of  the  work  on  the  analysis 
and  cancellation  of  singularities  in  Subsection  4.4.3. 

4.3  Polynomial  Entropies  for  the  Quadratic  Flux  System 

The  quadratic  flux  system  (2.5)  admits  entropy  functions  that  are  homoge¬ 
neous  polynomials  in  the  physical  variables  u  and  v  of  arbitrary  high  degrees. 
This  was  observed  in  [K]  for  the  symmetric  case  of  (2.5),  i.e.,  when  6  =  0.  We 
now  generalize  this  result  for  general  a  and  6.  We  remark  that  the  simple  func¬ 
tion  u 2  +  v2  is  a  strictly  convex  entropy  function  for  (2.5)  for  all  a  and  6.  This 
function  plays  a  special  role  in  obtaining  H~x  compactness  estimates  for  the 
dissipation  of  entropy  for  approximate  solution  sequences. 

Using  (4.1.2)  and  the  form  of  the  flux  matrix  in  (2.5),  we  obtain  the  entropy 
equation  for  (2.5)  as  follows: 


(4.3.1) 


((a  -  l)u  +  bv)r}uv  4-  (bti  +  v)(tjvv  -  tjuu )  =  0. 
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Propositon  4.3.1.  Given  any  a  and  b,  there  exists  an  infinite  sequence 
of  solutions  to  (4.3.1),  rjk  (k  =  1,2, 3,  •  -  •)  that  is  a  homogeneous  poly¬ 
nomial  in  u  and  v  of  degree  k. 

We  look  for  special  entropy  functions  of  the  form 

fc,,  x  « 

Vk  =  v  <p(a ),  a  = 

v 

Then  (4.2.1)  reduces  to  an  ordinary  differential  equation  for  4>(a).  It  is  then 
easy  to  check  that  for  each  positive  integer  k,  this  equation  admits  a  polynomial 
solution  in  a  of  degree  k.  Thus,  rjk  is  a  homogeneous  polynomial  in  u  and  v  of 
degree  k. 


4.4.  The  Riemann  Function  for  the  Entropy  Equation 


Next,  we  study  the  general  properties  of  entropy  in  the  (uq ,  iu2)-coordinates. 
In  particular,  we  are  interested  in  understanding  the  possible  singularities  of  rj 
and  the  construction  of  C 2  Goursat  entropies.  Goursat  entropies  axe  solutions 
of  the  entropy  equation  subjected  to  characteristic  boundary  conditions.  To 
this  end,  we  firs*  study  the  Riemann  function  11  for  the  hyperbolic  equation 
(4.1.3).  Recall  that  the  Riemann  function  contains  all  information  about  the 
general  solution  and  is  defined  as  the  solution  to  the  following  Goursat  problem 
(characteristic  boundary  value  problem)  with  special  boundary  data: 


U  =  K(wi,w2-,cr,T), 


n 


XJ)2 


+ 


W  2  —  W 1 


W2  —  W  i 


u>i 


=  o, 


(4.4.1) 


Huii  ^^(I.)  1 

- —  ds  >  , 

Urn  ^4  (£) 

~Z_y  ^  J  ' 


Consider  the  general  Goursat  problem  for  the  entropy  equation  (4.1.4) 


Vfw \W2  + 


M%) 


tu2  —  IV 1 

*?(G^2)  =  4>{W2), 


U>1 

w2  ~  W1 


~*lw\  — 


r)(wUK)  =  0(u>l), 


(4.4.2) 
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where  <f>  and  9  are  prescribed  functions,  and  f  and  k  fixed  constants.  It  is  well 
known  that  the  solution  rj  to  (4.4.2)  can  be  expressed  in  integral  form  in  terms 
of  TZ,  <j>,  and  9  (at  least  in  regions  where  T]  is  C2)  as  follows: 


/w a  /  -4i(^)  \ 

ft(u;i,u;2;f,s)  (  <t>'(s)  -  -j— ~  <j>(s )  J  ds 

+  J  n(wi,w2;y,K)  (o'(y)  +  dy' 


(4.4.3) 


Our  strategy  is  first  to  analyze  the  existence  and  singular  behavior  of  71 
and  then  make  use  of  (4.4.3)  to  show  how  to  choose  general  classes  of  data 
<(>  and  9  to  cancel  these  singularities.  The  existence  of  71  is  nontrivial  as  the 
coefficients  of  (4.4.1)  are  singular.  In  the  case  of  strictly  hyperbolic  systems,  the 
corresponding  7Z  and  77  equations  will  have  regular  coefficients  only  and  existence 
of  general  solutions  follows  from  standard  iteration  methods.  This  is  not  the  case 
here  and  is  one  of  the  key  difficulties  in  treating  nonstrictly  hyperbolic  systems. 
Our  analysis  on  7Z  and  77  will  be  applicable  to  more  general  situations  than 
(4.4.1)  -  (4.4.3).  To  state  our  results  in  their  strongest  form,  we  consider  the 
following  Goursat  problem  of  which  (4.4.1)  —  (4.4.3)  are  special  cases. 


(4.4.4) 


Vw,  It) 2 


+ 


B(wuw2 ) 

U72  —  Wi 


C{wx,w2) 

W2  —  Wi 


Vwi 


=  0, 


l(£,,w2)  =  <f>{w2), 


T](wUk)  =  9(wi). 


Here,  B  and  C  are  real  analytic  in  wx  and  w2  except  at  the  origin  (0, 0).  The 
corresponding  Riemann  function  satisfies 


71  =  71(wi,w2-,(t,t), 


7^w\  u>2  4" 


B{wi,w2)f>  C(wx,w2) 

"T 


71{w\,w2\ u;i,r)  =  exp 


"U>2 

W2  —  W 1  W  2  —  1^1 

**  B(wu») 


nwi  =  0, 


{-£ 


s  —  w  1 


dsj  , 


>  /  r  a  \ 

7l(wi ,  w2;  <t,w2)  =  exp<  -  /  - dy  >  . 

I  Jo  w2-y  J 


(4.4.5) 
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Theorem  4.4.1.  Consider  (4.4.5)  under  the  assumptions 

(1)  B  and  C  are  real  analytic  everywhere  except  at  (wi,w2)  = 

(0,0); 

(2)  There  exists  M  >  0  such  that  \B\oo,  |C|oo,  |(iui,u>2)V25(u>i,it>2)T|oo, 
and  \(wi,W2)V2C(wi,W2)t\oo  are  bounded  by  M  on  any  compact  set 

(J  =  1,2); 

(3)  k  =  0. 

Then  there  exists  a  unique  solution  'R(wi,w2‘,  £,0)  to  (4.4.5),  with  a  = 

£  and  r  =  k  =  0,  which  is  real  analytic  in  the  transformed  variables  w\, 
andw2  except  on  the  line  uq  +w 2  =  1,  werew 2  =  andwi  =  . 

The  proof  of  Theorem  4.4.1  relies  on  a  majorization  process.  A 
majorant  problem  corresponding  to  (4.4.5)  will  be  constructed  and  the 
existence  and  analyticity  of  'R,  follows  from  that  of  the  majorant. 

Lemma  4.4.2.  Under  the  assumptions  (1)  and  (2)  in  Theorem  4.4.1, 
there  exists  M1  >  0  such  that 

B{wx,w2)  ^  M ’ 

1  —  U?!  —  u72  1  —  uJi  —  w2  ’ 

C{wuw2)  M' 

l  —  W\  —w2  1  —  W\  —  w2 

Here,  <C  denotes  majorization  in  the  sense  of  real  analytic  power  series. 


Proof  of  Lemma  4.4.2.  We  will  prove  that 


(4.4.6) 


B(wl,w2)  M' 

1  —  Wi  —  W2  1  —  Wi  —  w2 


The  other  majorizing  inequality  involving  C  can  be  proved  in  a  similar  fashion. 


By  assumption  (1)  in  Theorem  4.4.1,  B{w\,w2)  is  real  analytic  everywhere 
except  at  the  point  (u>i,W2)  =  (1,0).  Thus,  away  from  this  point,  B  admits  a 
power  series  representation  of  the  form 

B(w\,w2)  =  ^2  cmnwimw2n. 

m,n>0 
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Moreover,  in  the  region  |uq  4-  w2\  <  1?  we  have 

1  _  I'm  +  n\ 

l-wx  -w2  jHE  \  n  ) 


wimW2n. 


m,n>  0 

Combining  the  two  power  series,  we  obtain 
B(w  i,w2) 

wx  —w  2 


5{W1,W2)  _  yr 

1  —  uJi  —  Wo  ■“ 


(-ran®!  w2 


(4.4.7) 


m,n>0 

=  ^  WiPW2g 
p,g>  0 


”E(mr 


^  W\mW2n 


0  <m<p  '  7 

0<n<ij 


The  formula  (4.4.7)  is  valid  in  the  region  |t/q  +  7Z72 1  <  1.  Using  (4.4.7),  we  trans¬ 
late  the  majorization  inequality  (4.4.6)  into  an  equivalent  inequality  involving 
the  coefficients  cm,n  in  the  power  series  expansion  of  B, 

B{wx,w2)  M' 

1  —  Wi  —  w2  1  —  Wi  —  w2 


\  y  / 


p>0 

0<n<g 


(4.4.8) 


5 


p>  o 
0<n<q 


Our  task  is  to  verify  the  last  inequality  in  (4.4.8).  To  this  end,  we  consider 
the  auxiliary  function  ,  where  m  is  an  integer  and  g(z)  =  Ylk>i  9kZk  is 

an  analytic  function  of  2  in  the  disk  |z|  <  1.  We  have  the  following  formulae: 

9(z) 


(1-z) 


m+1 


(4.4.9) 


=e^‘e(T>“ 

k>\  n>0  X  7 

=E2'  e  (m:v 


g>l  0<n<g 


We  now  make  a  useful  choice  of  g(z)  that  is  related  to  B.  First,  we  rewrite 
B(w\,w2)  in  the  following  way: 

B(wi,w2)  -  ^ gi(w2)wi\ 
i>  o 

(4.4.10)  gi(w2 )  =  ^  clkw2k. 

k>0 
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We  will  take  g  to  be  gi  in  (4.4.9).  We  claim  that,  on  |itf2|  <  1, 


(4.4.11) 


gP-m(w2)  M 

(1  -  uJ2)m+1  (1  -ttJ2)m+1  ’ 


where  M  is  a  constant  independent  of  p  and  m.  The  proof  of  the  claim  de¬ 
pends  on  certain  decay  estimates  of  the  coefficients  cmn.  Consider  B(poet6 ,  z)  = 
Yli> o  9i(z)Poette  where  0  <  po  <  1  and  z  belongs  to  the  unit  disk  on  the  complex 
plane.  B(poetff,z)  is  analytic  for  all  \z\  <  1.  Moreover,  we  have 


gi(z)  =cp0l  /  B(pQe,6,z)e  ,lB  dO, 
Jo 


\gi(z)\  <  2ircp0l\\B\\oo. 

Here  c  is  a  universal  constant.  Now  let  pQ  — >  1“,  we  obtain 


(4.4.12) 


lsi(*)l  <  2«||B|U. 


Likewise,  we  have 

9i'W  =cp0-'  f  d!,B(p0e‘\z)e-iU  d6. 

Jo 

Using  the  bound  \\z2dlB(poe'e ,  ^)lloo  <  c',  where  c'  is  independent  of  po  and  6 , 
we  obtain  in  a  similar  fashion, 

(4.4.13)  ||*V(*)lloo  <  27rc", 

where  c"  is  a  constant  independent  of  /. 

Combining  (4.4.12)  and  (4.4.13),  it  is  not  difficult  to  obtain  a  decay  estimate 
for  the  coefficients  of  B 

(4.4.14)  M  <  M 
Now  (4.4.11)  is  equivalent  to 
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The  right  side  satisfies 


M 


' m  +  q N 


V  9  /  (9  -  "  +  l)2 

'  ^  '  n>9 

Efm  +  n\ 

.A  «  > 


-p  — m,g  — n  | 


0<n<? 


This  proves  (4.4.11). 


We  now  show  that  (4.4.11)  and  (4.4.14)  implies  this  lemma.  Using  (4.4.10) 
and  (4.4.9),  we  see  that  (4.4.11)  is  equivalent  to 


(4.4.15)  £ffij'  £  (m+n)cP-m.,-n^  £ 


M 


g>0  0  <  m  <p 

0<n<<j 


0<m<p 


(1  -u72)m+1' 


The  right  side  of  (4.4.15)  can  be  summed  to  give 

M  v-''  /p  +  9  +  1N 


(4.4.16) 


£ 


9  +  1 


w2 


(4.4.16)  is  now  easily  checked  to  imply  (4.4.8)  by  virtue  of  (4.4.14). 

This  completes  the  proof  of  Lemma  4.4.2. 

Using  Lemma  4.4.2,  we  obtain  a  majorant  problem  of  (4.4.5). 

Theorem  4.4.3.  Suppose  G  =  G{w  1,102)  is  real  analytic  in  a  domain 
Q  and  that  Q  is  the  solution  to  the  following  majorant  problem: 

G  =  G{w  i,u>2), 

M' 

Gwi  u>2  —  Z - - - ~(Gwi  +  Gw 2)  —  0, 

I  —  W  i  —  W  2 


G(wu0)  = 


M' 


1  —  ’ 

1  —  Uq 


M' 

G(0,W2)  =  - — • 

1  —  W2 


(4.4.17) 
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Then,  if  7Z  is  the  solution  to  (4.4.5)  with  r  =  k  =  0,  we  have,  on  f l, 

TZ  ■<  Q  in  the  sense  of  power  series  in  w\  and  W2- 

We  remark  that  (4.4.17)  is  the  celebrated  Euler-Poisson-Darboux  equation. 
We  summarize  some  of  its  most  fundamental  properties  in  the  next  lemma.  We 
also  note  that  although  the  analyticity  of  Q  (see  Lemma  4.4.4  below)  implies  that 
of  7 Z,  71  is  not  everywhere  analytic  on  the  domain  of  interest  to  our  problem. 
It  turns  out  that  the  line  w\  =  W2  (and  hence  the  umbilic  point)  lies  on  the 
boundary  of  the  region  of  convergence  of  the  power  series  for  71  and  Q.  On  this 
part  of  the  boundary,  7Z  and  its  derivatives  are  singular.  We  characterize  these 
singularities  in  Lemma  4.4.4  and  Lemma  4.5.1. 

Lemma  4.4.4. 

(1)  In  the  variables  Wi  =  1  —  w\  and  W2,  (4.4.17)  is  an  Euler- 

Poisson-Darboux  equation  with  special  Goursat  data: 


Gu)\W2  _  -  (^Wl  Gw?)  —  0? 

W  2  —  10 1 

M' 

a  ^i,o  =  — , 


(4.4.18) 


Q{0,w2)  = 


1  —  w2 


(2)  The  Riemann  function  of  the  first  equation  in  (4.4.18)  is 


n(ivi,w2;$,K)  =  7~. - =-tt77 F{M', M',  l;p), 


(4.4.19) 


(u;i  -  w2)2M' 

(wi  —  k)M'(£  —  W2)M' 
(£  -  W\){W2  -  «) 

(wi  -  k)(z  -  w2y 


Here,  F  is  the  Hypergeometric  function  defined  by  the  power  series 

rv  a  ^  ,  ,  f  fl(a  +  l)-(a+m-ip  +  l)--(^  +  m-  l)>B 

F(a,/37;p)  =  l+^  - 1.2...m.T(7+1)...(T  +  m_1) - "  ' 


(4.4.20) 


\P\  <  1- 
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(3)  The  solution  to  (4.4.18)  is  given  by 


/Wi 

'H(wi,W2;£,0)M'(M'  -  1)£-2  d£ 

(4.4.21)  +  /  H(wi,W2\  l,y)M'(l  —  M')(l  —  y)~2  dy. 

Jo 


The  proof  of  Theorem  4.4.1  follows  from  Theorem  4.4.3  and  Lemma  4.4.4. 
4.5.  C 2  Goursat  Entropies 

In  this  subsection,  we  characterize  the  singularities  of  entropies  and  we  show 
how  to  construct  Goursat  entropies  that  are  regular.  We  continue  to  work  on 
the  general  level  of  equations  (4.4.4)  and  (4.4.5). 

Similar  to  (4.4.3),  we  have  an  integral  representation  of  fj  in  terms  of  71. 

v(w1,W2)  =  0(Z)TI(wi,W2;Z,k)+  [  H(wl,w2;t,s)(p'(s)  +  ~^,S}(l){s))ds 

Jk  s-Z 

(4.5.1)  +  f  H{wl,w2;y,K)(d'(y)  +  ^^-6{y))dy. 

Js  K-y 


Using  (4.5.1)  with  k  =  0,  it  is  clear  that  whatever  singularities  fj  might 
have  are  inherited  from  those  of  71.  Since  1Z  C  Q,  the  singularities  of  71  are 
dominated  by  those  of  Q .  We  make  this  precise  with  the  aid  of  the  next  lemma. 


Lemma  4.5.1.  Suppose  that  Fi  and  F2  are  real  analytic  on  an  open 
domain  Cl  C  R2.  Assume  further  that 

(1)  Fi  <  F2  on  Cl, 

(2)  F2  —  ^F2  where  F2  is  analytic  on  Cl  and  S  is  analytic  on  Cl. 
Then  we  have 


Fi  =  ~FU  where  Fx  =  SF, 
o 


is  analytic  on  Cl. 
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We  have  characterized  the  singularities  in  Q  and  its  derivatives.  By  repeated 
applications  of  Lemma  4.3.1,  we  can  show  that  the  singularities  of  7Z  and  its 
derivatives  are  dominated  by  those  in  Q  and  corresponding  derivatives.  Thus,  we 
have  established  a  kind  of  comparison  principle  for  the  singularities  of  solutions 
to  second  order  hyperbolic  equations  of  Euler-Poisson-Darboux  type  One  of 
our  main  objectives  is  to  construct  regular  solutions  to  (4.4.4).  It  is  now  clear 
that  this  can  be  achieved  if  we  can  prescribe  Goursat  data  that  could  can  eel  the 
type  of  singularities  in  Q  since  they  dominate  the  singularities  in  71  which  is  the 
Riemann  function  for  (4.4.4). 


We  now  make  use  of  what  have  been  developed  so  far  to  construct  general 
classes  of  entropies  that  are  C 2  in  w\  and  W2 .  The  domain  of  interest  for  our 
problem  is  {io i  <  0  <  W2}.  It  is  enough  to  focus  our  attentions  on  a  compact 
set  containing  the  umbilic  point  <  uq  <  0  <  u>2  <  w£ }  where  >  0  and 
rej”  <  0  are  fixed  constants.  Also  fix  a  constant  8  >  0  satisfying  |£|  <C  w £ ,  — irf. 


Without  loss  of  generality,  we  consider  (4.4.4)  with  the  special  choice  of 
Goursar  data  <f>(w2)  =  0: 

^  ,  B(wuw2)_  ,  C(wuw2)_  _n 

w?  1  TJw?  r  _  _  Tjw I  —  U, 


(4.5.2) 


W  2  —  Wi 
n(€,w  2)  =  0, 

fi(w  1,«)  =  B(wi). 


W2  -  w  1 


There  is  no  loss  of  generalities  since  an  analysis  similar  to  what  is  shown 
below  will  allow  us  to  construct  C 2  entropies  with  Goursat  data  9  —  0  and  <f> 
arbitrary  (up  to  a  finite  number  of  conditions,  see  Theorem  4.5.2).  The  general 
case  (as  in  (4.4.4))  then  follows  from  a  simple  linear  superposition  of  these  two 
special  cases. 

Using  the  fact  that  TZ  C  Q  and  combining  Lemma  4.4.4  and  Lemma  4.5.1, 
we  obtain  a  decomposition  of  It  into  a  product  of  a  singular  and  a  regular 
part.  The  regular  part  is  analytic  everywhere.  A  similar  decomposition  exists 
for  Q  using  Lo’inma  4.4.4.  Moreover,  from  Lemma  4.5.1  again,  it  is  possible 
to  set  the  singular  parts  to  be  identical  for  both  It  and  Q.  Then  using  this 
decomposition  with  (4.5.1),  it  is  possible  to  obtain  conditions  on  9  that  guarantee 


the  cancellation  of  singularities.  We  now  state  our  main  result  in  this  section. 


Theorem  4.5.2.  Given  a  positive  integer  m  and  given  a  compact  set 
K  containing  (0,0)  in  R2.  Consider  the  Goursat  problem  (4.5.2)  with 
k  =  0.  There  exists  a  subspace  M  C  X2[u>f,  — <5]  satisfying  codimM  < 
oo  such  that  if 

(1)  deMnCm, 

(2)  6(w\)  =  0,  for  all  w\  >  —6, 

then  there  exists  solution  fj  £  Cm([ioj” ,  0]  x  [0,  wf  ])  to  (4.5.2)  satisfying 
max  sup  \\wi\~k\w2\~ld^dlW2fi\  <C. 

0<k+l<m  (WliU,2)£/i 

Recall  that  (4.4.1)  —  (4  ‘.3)  axe  special  cases  of  (4.4.4), (4.4.5)  and  (4.5.1). 
We  obtain  immediately  existence  of  many  Cm  Goursat  entropies  for  the  qua¬ 
dratic  flux  sj^em  (2.5). 

Corollary  4.5.3.  Under  the  same  assumptions  on  the  Goursat  data 
as  in  Thee  <  ?m  4.5.2,  there  exists  a  solution  n  €  Cm([ie(~ ,  0]  x  [0,u>^"]) 
satisfying 

max  sup  \\wi\-k\w2\~ld^ldlWjT)\  <C. 

°<k+l<m(WllW2)€K 

We  now  turn  to  investigate  the  regularity  of  entropies  constructed  in  The¬ 
orem  4.5.2  and  Corollary  4.5.3  as  function  of  the  state  variables  u  and  v.  This 
regularity  is  needed  in  our  use  of  compensated  compactness  to  prove  the  con¬ 
vergence  of  approximate  solutions  to  (2.1). 

In  many  strictly  hyperbolic  systems,  such  as  the  system  of  elasticity,  the 
map  J  :  ( u,v )  — »  (101,102)  is  C2.  Therefore.  C2  regularity  of  entropies  in  the 
(u,  v)-coordinates  is  a  direct  consequence  of  Theorem  4.5.2.  For  nonstrictly 
hyperbolic  systems,  the  coincidence  of  eigenvalues  Aj  and  A2  usually  means  that 
the  geometry  of  the  wave  curves  are  very  singular  at  the  umbilic  point.  Thus, 
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J  is  usually  not  C 2  and  additional  work  is  necessary  to  get  C 2  regularity  for 

Tj  =  T](U,V). 

We  now  restrict  ourselves  to  region  IV  of  the  quadrat’c  flux  system  (2.5), 
i.e.,  where  a  and  b  satisfy  A  =  —  3264  +  62[27+36(a— 2)  —  4(a  —  2)2]+4(a— 2)3  >  0. 
Using  Proposition  3.1,  we  have 

Proposition  4.5.4.  Consider  the  quadratic  £ux  system  (2.5)  in  region 
IV,  i.e.,  A  >  0.  Near  the  u mbilic  point  (u,  v)  =  (0,0),  the  derivatives 
of  the  Riemann  invariants  satisfy  the  following  estimates: 

wt  =  0(1),  i  =  l,2, 

71)  ' 

(4.5.3)  sra>i  =  O(-p^),  1  <  m  +  n  <  2. 

Next,  we  recover  from  the  proof  of  Theorem  4.5.3,  certain  fine  estimates 
of  the  derivatives  rj  in  Wi  and  W2-  As  we  shall  see,  these  derivatives  vanishe  to 
certain  orders  as  they  approach  the  umbilic  point. 

Proposition  4.5.5.  Consider  (2.5)  in  region  TV.  Given  any  nonnega¬ 
tive  integers  m.  n ,  k,  l  with  m  >  k,n  >  l,  there  exists  entropy  functions 
constructed  using  the  method  in  Theorem  4.5.3  and  which  satisfy  the 
following  estimates  near  the  umbilic  point: 


Now  using  the  chain  rule  and  combining  the  estimates  in  Propositions  4.5.4 
and  4.5.5,  we  obtain 

Theorem  4.5.6.  Consider  (2.5)  with  a  and  b  satisfying  A  =  —32b4  + 
b2(27  +  36(a  —  2)  —  4(a  —  2)2)  +  4(a  —  2)3  >  0.  Let  J  :  (u,  v)  — ►  (uq,  ic2) 
denotes  the  map  from  the  state  space  to  the  Riemann  invariants  plane. 
Suppose  that  r?  is  the  solution  to  (4.4.2)  as  constructed  in  Corollary 
4.5.3.  Then  we  have 


T)oJe  C2(R2). 
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5.  Approximate  Solutions  and  Young  Measures 
5.1.  Parabolic  Approximations 

We  consider  parabolic  approximations  to  the  general  system  (1.1)  by  adding 
artificial  viscosities.  Recall  that  (1.1)  together  with  initial  data  Uo(x)  takes  the 
form 


dtU  +  dxF{U )  =  0,  t  >  0,2  €  R, 
(5.1.1)  U(x,Q)  =  Uo{x). 


Consider  a  sequence  of  parabolic  approximate  solutions  Ue  governed  by  the 
associated  parabolic  system 

dtU(  +dxF(Ut)  =  edx(DdxU(),  t  >  0,x  €  R, 

(5.1.2)  U*(x,  0)  =  Uo(x), 

where  D  >  0  is  a  non-negative  matrix  (viscosity  matrix)  and  e  >  0  measures  the 
amount  of  artificial  viscosities  in  (5.1.2). 

We  are  concerned  with  the  convergence  of  Uf  to  a  weak  solution  U  of  (5.1.1) 
as  e  — ♦  0+.  A  compactness  framework  theorem  (Theorem  6.3)  will  be  established 
in  Section  6  to  attain  this  goal.  It  will  then  be  applied  to  the  quadratic  flux 
system  (2.5).  This  is  made  possible  by  a  reduction  of  Young  measure  analysis 
(Theorem  6.2)  using  large  classes  of  regular  entropy  functions  (Theorem  4.5.2 
and  Theorem  4.5.3)  constructed  in  Section  4.  In  this  section,  we  dispense  with 
certain  preliminaries  concerning  existence  and  L°°  apriori  estimates  of  solutions 
to  (5.1.2).  For  simplicity,  we  take  D  =  7,  the  identity  matrix. 

Proposition  5.1.1.  Let  ci  and  c2  be  constants.  Let  0ClC2  =  {(u,  v)  |  ci  < 
uq  <  0  <  w?  <  C2}.  Suppose  that  the  Riemann  invariants  w \  and  w ? 
are  quasi-convex  on  dQClC2.  That  is,  on  dQClC2> 


r[  •  V2w2  ■  ri  >  0, 

-  •  V2uq  •  r2  >  0. 
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Then  that  QClC2  is  an  invariant  region  for  (5.1.2),  i.e.,  Uq(x)  €  flClc2, 
for  all  x  €  R,  implies  Ue(x,t )  6  fiqcj)  for  all  x  €  R,  f  >  0.  Moreover, 
if  (5.1.2)  admits  a  family  of  such  invariant  regions  which  spans  R2, 
then  we  obtain  an  apriori  L°°  bound  for  solutions  Ue  to  (5.1.2). 

Proposition  5.1.1  is  a  straightforward  consequence  of  the  invariant  region 
theorem  (see  [CCS]).  Granting  such  an  L°°  apriori  estimate,  a  local  in  t  solution 
to  (5.1.2)  obtained  by  standard  iteration  arguments  can  be  extended  globally  in 
t  >  0. 

5.2.  Finite  Difference  Approximations 

As  is  known,  it  is  very  complicated  to  directly  construct  the  Riemann  so¬ 
lutions  for  nonstrictly  hyperbolic  systems  of  conservation  laws  (cf.  [IMPT,  IT, 
SS2,  SSMP]).  Fortunately,  we  can  obtain  the  global  entropy  solutions  of  Riemann 
problems  in  any  bounded  domain  in  (x,t)-plane  by  proving  the  convergence  of 
the  viscous  approximations  and  noting  the  finiteness  of  propagation  speeds  of 
the  entropy  solutions. 

Now  we  construct  Lax-Friedrichs  approximations  [La3]  and  Godunov  ap¬ 
proximations  [Go]  Ul(x,t )  =  ( ui(x,t),vl(x,t )). 

On  the  upper-left  plane  t  >  0,  we  have  the  grid 

t  =  nh  ,  x=jC]  j  =  0, ±1,±2, •••  ,  n  =  0, 1,2, , 

where  n  >  0  and  j  are  integers  and  the  positive  constants  h  and  i  are  the  time 
step  length  and  the  space  step  length,  respectively,  that  satisfy  the  inequality 

£ 

(5.2.1)  max(  sup  ]A,(£/*)|)  <  —  <  M  , 

*=  1  j2  —  oo<x<oo  tl 

0  <t<T 

for  any  given  T  >  0.  We  shall  prove  that  Ue(x,  t)  are  uniformly  bounded  that 
(5.2.1)  always  holds. 


5.2.1.  Lax-Friedrichs  Approximations 


30 


For  integers  n  >  1,  we  set 

Jn  =  {j  •  j  integers,  n  +  j  =  even}  . 


On  the  rectangle  {(x,t)  :  (j  —  l)i  <  x  <  ( j  +  1)^,  0  <  t  <  h,  j  odd},  we 
define  Ut{x,  t)  as  the  solutions  of  Riemann  problems 


r  (i-i)  , 


< 


U\t=o  = 


I  U‘((j  -  1)£)  , 

1.^(0+ m, 


x  <  j£  , 
x  >  jl , 


where  Uq(x)  =  Uo(x)Xe(x),  and 


and  define 


Suppose  that  (ue,  ve)  have  been  defined  for  t  <  nh.  Then,  on  the  rectangle 
{(x,  t)  :  j£  <  x  <  ( j  +  2)^,  nh  <  t  <  (n  +  l)h,  j  £  J„},  we  define  Ul(x,  t )  as  the 
solutions  of  Riemann  problems 


and  define 


(5.2.2) 


'(1-1), 

-  ...  x<(j  +  ik, 

,  {Ufa,  x>(j  +  l)t, 

1  rU+l)t 

U"+'  =  ±.  U‘(x,(n  +  l)h-0)dx  . 
"  J(j- IK 


In  this  manner  we  obtain  the  Lax-Friedrichs  approximations  on  the  upper-half 
plane,  and  the  difference  scheme  (5.2.2)  is  just  the  Lax-Friedrichs  scheme  [La3]. 


5.2.2.  Godunov  Approximations 
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Similarly,  on  the  rectangle  {(x,t)  :  j£  <  x  <  (j  +  1)£,  0  <  t  <  h},  we  define 
Ut(x,t)  as  the  solutions  of  Riemann  problems 


'(1-1), 

•  r  usui) .  *  <  a + |k , 

.  1,-0  \v‘(U  +  W,  *>0'  +  |K. 


where  Uq(x)  =  Uo(x)Xt(x),  and 


and  define 


Xt(x) 


j1- 

lo. 


v-=\L 


(i+l/2)< 


(j-1/2)* 


otherwise, 

[/*  (x,  h  -  Q)dx  . 


Suppose  that  Ue  have  been  defined  for  t  <  nh.  Then,  on  the  rectangle 
{(x,t)  :  j£  <  x  <  ( j  +  1)£,  nh  <  t  <  (n  +  l)/i),  we  define  Ul(x,t)  as  the  solutions 
of  Riemann  problems 


■(in . 

ii-.»  = 


u?  ,  X  <  (j  +  \)(  , 

Ufr,,  *  >  u  +  |K . 


and  define 

1  fO+l/2)* 

(5.2.3)  C/.n+1  =  7  /  tf*(x,  (n  +  l)fe  -  0)dx  . 

This  completes  the  construction  of  the  Godunov  approximations. 


Remark  1.  The  approximate  solutions  U((x,  t )  constructed  above  have  the  same 
local  structure  as  the  random  choice  approximations  of  Glimm  [G2]. 


Remark  2.  If  the  Riemann  solutions  have  convex  invariant  regions  in  U- space,  it 
is  easy  to  check  that  the  Lax-Friedrichs  and  Godunov  approximations  have  the 
same  convex  invariant  regions  in  U- space  from  their  construction  and  Jensen’s 
inequality.  This  enables  us  to  obtain  L°°  uniform  bounds  for  these  approxima¬ 
tions. 
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5.3.  Young  Measures  and  Compensated  Compactness 

The  Young  measure  representation  (see  [Ta])  for  the  sequence  of  bounded 
functions  in  an  appropriate  space  is  am  efficient  tool  for  studying  the  limit  be¬ 
havior  of  the  approximate  solutions  of  nonlinear  problems,  especially  for  conser¬ 
vation  laws  because  of  the  lack  of  regularity  of  the  limit  problems.  By  combining 
the  Young  measure  representation  with  compensated  compactness  [Ta.Mu],  one 
can  transfer  the  singular  limit  problems  to  the  problems  of  solving  some  func¬ 
tional  equations  for  the  corresponding  Young  measures,  that  is,  to  studying  the 
structure  of  the  Young  measures  satisfying  the  functional  equations.  If  one  can 
solve  these  functional  equations  to  clarify  the  structure  of  the  Young  measures, 
the  limit  behavior  of  the  sequence  can  be  well  understood.  Therefore,  the  essen¬ 
tial  difficulty  is  how  to  solve  these  functional  equations  for  the  Young  measures. 
This  difficulty  is  overcome  for  some  important  systems  in  conservation  laws  (cf. 
[Di,  Chi,  Ch2,  DCL,  K,  Mo,  Se]).  In  this  section  we  review  some  results  on  the 
Young  measures  and  the  compensated  compactness  for  subsequent  development 
to  solve  general  nonstrictly  hyperbolic  systems. 

Theorem  5.3.1([Taj).  Suppose  that  Uf  :  R+  — ►  Rn  is  a  sequence  of 
bounded  measurable  functions 

(5.3.1)  U'(x,t)  €  K,  a.e. 

for  a  bounded  set  K  in  Rn.  Then  there  exists  a  subsequence  (still 
labeled  Ut)  and  a  family  of  Young  measures 

vx,t{ A)  €  Prob.(Rn), 


such  that 

(1)  (i)  for  any  continuous  function  g, 

w*  -  lim  g(Ut)  =<  i'i,t(^)>ff(^)  >=  /  sr(A)<fi/Xi,(A); 

JRn 

(2)  (ii)  Uf(x,t)  — ►  U(x,t )  strongly  if  and  only  if  vZtt  is  a  Dirac  mass 


vx,t  =  &U(z,t), 
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for  almost  all  (x,  t); 

(3)  (iii)  uXtt  satisfies 

<  ux,tiT}\  >  <vx,t,qi> 

<vz,t,q2> 


(5.3.2) 


<  vx 


m  qi 

T}2  92 


provided  that 

(5.3.3)  T]i(Ut)t  +  qi{Ut) x  compact  in 

for  continuous  function  pairs  (rji,qi),i  =  1,2. 


This  theorem  ensures  the  existence  of  the  Young  measures  uniquely  deter¬ 
mined  by  the  sequence  of  bounded  functions.  The  second  result  indicates  that 
the  strong  convergence  of  the  sequence  is  equivalent  to  the  one  point  structure 
of  the  support  of  the  Young  measures.  Furthermore,  the  boundedness  of  Uf 
(5.3.1)  automatically  ensures  that 

for  continuous  function  pairs  ( r}i,qi),i  =  1,2.  This  theorem  indicates  that  an 
extra  condition  of  weak  composite  compactness  for  rji(Ut)t  +  qi{Ut)x  can  give 
us  very  useful  information  for  the  Young  measures.  These  theorems  provide  a 
framework  by  which  one  can  prove  strong  convergence  of  the  sequence  Ut(x,t) 
satisfying  (5.3.1)  and  (5.3.2)  by  deducing 

from  the  functional  equations  (5.3.3)  for  some  continuous  function  pairs. 

The  following  compactness  embedding  theorems  are  useful  for  obtaining  the 
condition  (5.3.2)  for  conservation  laws. 

Theorem  5.3.2.  Let  l<q<p<r<oo.  Then 


( compact  set  of  H  ( bounded  set  of  Wj0J.’r)  C  ( compact  set  oflVjo1c’p). 
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The  proof  of  Theorem  5.3.2  can  be  found  in  [DCL,  Ch2]. 

Theorem  5.3.3([Mu]).  The  embedding  of  the  positive  cone  of  W~1,p  in  W~1,q 
is  completely  continuous  for  q  <  p. 

Theorem  5.3.2  indicates  that  compactness  in  coupled  with  bound¬ 
edness  in  yields  compactness  in  Theorem  5.3.3  says  that  the 

uniformly  lower  (or  upper)  bound  in  the  dual  sense  in  W~l,p  of  the  sequence  in 
W~l,q  leads  to  compactness  in  W~1,q,q  <  p. 

5.4  The  Dissipation  Measures  77(1/*)*  +  q(Ut)x 
5.4.1  Parabolic  Approximate  Solution  Sequences 


Consider  a  sequence  of  viscosity  approximate  solutions  {U(}e> 0  to  (5.1.2). 
Now  suppose  that  77,  is  a  C 2  strictly  convex  entropy  for  (5.1.1),  and  that  Uo(x ) 
tends  to  a  constant  state  U  as  |x|  — >  00  and  Uo  —  U  6  L2  fl  L°°.  In  the  case  of 
(2.5)  we  may  choose  77,  to  be  u2  +  v2.  Multiplying  (5.1.2)  by  Vrj^(Uf)  —  Vt7*(I7), 
a  standard  integration  by  parts  argument  gives  the  estimate 


\U(X\2  dx  dt  <  C, 


where  C  depends  on  Uq. 


Consider  any  C 2  entropy-entropy  flux  pair  (77,5)  for  the  system  (5.1.1).  In 
the  case  of  (2.5),  such  pairs  were  constructed  in  Section  4.  Multiplying  (5.1.2) 
by  V77,  and  integrating  by  parts,  we  get,  after  using  the  L°°  bound  on  U %  the 
boundedness  of  V2t]  on  compact  sets  and  a  standard  application  of  Theorem 
5.3.3  [Mu],  a  weak  compactness  estimate  for  the  dissipation  measure 

77 {U()t  -I-  q(Ue) x  €  compact  set  of  H^0\. 


By  the  L°°  bound  of  U(,  as  e  — ►  0+,  U(  converges  in  weak  star  topology  in 
L°° .  Then,  by  Theorem  5.3.1  [Ta],  there  exists  a  family  of  probability  measures 
u(z,t)i  the  Young  measures,  that  describes  this  weak  convergence  in  the  following 
way.  For  all  continuous  function  g ,  we  have 

w  *  —  limi^ Q+g(Uf(x,t))  = 


J  gW  <fyz,t)(A)=  <  >  • 
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Let  ( T],q )  and  (77,  q)  be  two  pairs  of  C 2  entropy-entropy  flux.  They  satisfy 
compactness  conditions  in  H~l  as  stated  above.  Then,  by  Theorem  5.3.1  again, 
we  get  the  commutation  relation: 

(C)  <  v,r\q-  fjq  >=<  v,  r)  ><  u,  q  >  -  <  v,  7?  ><  u,  q  >  . 

Here  we  have  dropped  the  subscript  (x,t)  on  V(z,t)- 


5.4.2  Finite  Difference  Approximate  Solution  Sequences 

Consider  approximate  sequences  {Z7*}*>o  generated  by  either  the  Godunov 
or  the  Lax-Friedrichs  schemes.  Using  Green’s  formula,  we  obtain  that,  for  any 

4>  €  Co(IIt), 


(541)  II 


where 


0  <t<T~mh 


(v(Ue)<f>t  +  q(Ut)(f>x)dx  dt  =  M(4>)  +  L{4>)  +  S(<f>)  , 


/oo  1 00 

4>(x,T)T)(U\x,T))dx  -  /  (f>(x,  Q)r)(Ue(x,  0))dx  , 

-00  J  —  CO 


fT 

(5.4.3)  S(<f>)  =  /  -  [ q ]}  dt  , 

Jo 


(5.4.4) 


For  the  Lax-Friedrichs  approximations,  n  +  j  —  even, 

r(j+i)* 


(5.4.5) 


=  /  (; n(uLn)-v(u‘n))dx  , 


'u- d< 


W)  =  E  /  (v(ut-)-Ti(utj’'))(4>-<f>';)dx 

It 


For  the  Godunov  approximations, 


(5.4.6) 


r(j+l/2)t 

W)  =  E*"/  (q(ub)-Ti(u'"))dx , 

“  70-1/2)' 

y(j+l/2)* 

w)  =  E  /  (^- )  -  ^n))  (*  -  *")  ^  . 

7^  Aj-1/2K 
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where  Uin  =  U^x^nh  —  0),  =  <j>(j£,nh),  the  summation  in  is  taken  over 

all  shock  waves  s  in  U(  at  fixed  time  t,  and  a  is  the  propagating  speed  of  the 
shock  wave.  If  s  =  (:r(£),£),  then  [77]  and  [5]  denote  the  jump  of  Tj(Ut(x,t))  and 
q(Ue(x,t ))  across  s  from  left  to  right  respectively,  namely, 

(  M  =  r]{U\x{t)  +  0,i)}  -  j]{Ue(x(t)  -  0,#)}  , 

\  [q]  =  q{ue{x(t )  +  0,  t)  -  q{ul(x(t)  -  0,  t)} . 


To  arrive  at  (4.4),  we  first  notice  that  U  have  compact  support  in  the  region 
nT  and,  therefore,  we  may  substitute  tj  =  77*,  q  =  <7*  and  <f>  =  1  in  the  equality 
(5.4.1).  Then 


(5.4.7) 

while 


E  /  MV-)  -  v*(u?))dx  +  /  E^M  -  Iff*]} 
=  J  ri.(Ue(x,fi))dx  -  J  q*(Ue(x,T))  dx 
<  f  n  vo(x))  dx  <  C  , 

J  —  OC 


dt 


_  tU+ 

E  /  MUl)  -  r,m(U*))dx 

hn  A;- Dt 

(5.4.S) 

=  E  /  dx  (l  -  9){U2  -  U")T  V2Tjm(U"  +  9{U1  -  U?))(UZ  -  U”)d9  . 
7 n  hi-' )«  A) 


Notice  that  the  entropy  inequality  crfr/*]  —  [g„]  >  0  is  satisfied  across  the 
shock  waves  and  tj,  is  a  convex  entropy.  We  obtain  from  (5.4.7)-(5.4.8)  that 

rT 

(5.4.9)  /  E -  \qA) dt  <  c  , 

Jo 

and 

(5.4.10) 

E  /  dx  (l  -  9){U2  -  +  9(U2  -  UJ))(U2  -  U”)d9  <  C  . 

j,n  Jo 
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In  particular,  since  rT  •  V2r;,  •  r  >  Co(r,r),  Co  >  0  constant,  we  obtain 


(5.4.11) 


2  -  Uf\2dx  <  C  . 


For  any  bounded  set  ft  C  IIt  and  weak  entropy  pair  (77,  q),  we  derive  from 
(5.4.2)-(5.4.5),  and  (5.4.9)-(5.4.11)  that 

\M(<j>)\  <  cmCom , 

\s(<i>)\  <  M\Co  [  E  _  MM* 

Jo 

<C\\4>\\CoJ^  5^{<r[r/,]  -  [5*]}  dt 

_  rU+ IK 

i^wi  =  iE^n  / 

<iwic.D/  (<K«)-i(p;))<fci 

77  Ai-i>* 


j,n 

_  /•(>+!)<  fl 

<wiic„E  /  /  (i-w--^)Tvv^?+^-^,))(^--^") 

Jo 

-  ^IMIIc0(n)  • 


d<? 


Hence 

\(M  +  Lx  +  S)(<j>)\  <  C\\<f>\\Co  ; 


that  is, 

WM  +  Li+SWc.  <C  . 


Using  the  embedding  theorem  Cg(ft)  W  1,70 (ft),  1  <  qo  < 


we 


have 


(5.4.12) 


M  +  L\  +  S  compact  in  IF  1,70 (ft)  . 
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Furthermore,  for  any  <f>  €  C'o'(O),  |  <  a  <  1,  we  have  from  (4.31)  that 
_  c(j+iR 

\LM) I  <  E  /  Wx,nfc)-^”|  |,(tr;)-,(t/;)|  dx 

hi-iv 

_  _  Hi+iK 

<  nwic?  DE  /  W-)  -  ^;)i2^  )1/2 

r(j+l)< 

<  ||v,||t_r-1/2iMU(E  /  |K  -  t/;i2<ix  )■/* 

j,„  Ai-D< 

<  cr-'/2||^||c;(a) 

<C£a  1/2|Wlliv„'  »(n)  ’  P>  ’ 

that  is, 

(5.4.13) 

l|i2lllv-..,.(n)<«“-1/2-.0,  («-.0),  1  <»<■--”- 

TX  X  i”  Ot  TX  X 

We  obtain  from  (5.4.12)  and  (5.4.13)  that 

M  +  L  +  5  compact  in  W~1,7o(fl), 
with  the  aid  of  Theorem  5.3.2  and  Theorem  5.3.3.  Therefore, 

(5.4.14)  v(Ul)t  +  q(Ul)x  €  compact  set  of  H j“*. 

Similar  to  Section  5.4.1,  we  can  obtain  the  commutation  relation  (C)  from 

(5.4.14) . 

6.  Compactness  Framework  for  Approximate  Solutions 

Now  we  establish  a  structure  framework  for  the  Young  measures,  which  are 
uniquely  determined  by  the  approximate  solutions  to  the  system  (1.1)  with  an 
isolated  umbilic  point  P  =  ( wi,W2 ): 

Ai(P)  =  A2(P), 

in  the  Riemann  coordinates.  Then  we  conclude  a  corresponding  compactness 
framework  for  approximate  solution  sequences  to  the  system  (1.1). 

First,  we  list  some  basic  assumptions  on  the  structure  of  (1.1)  for  the  frame¬ 
work  theorems  that  we  will  state. 
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Hypothesis  on  the  Hyperbolic  System  (1.1)  for  Theorem  6.1. 

/  TT-.  \  A'U>,  Aj(wUW2) 

{  )  A2-Aj  W2-W,  ’ 

where  i  ^  j,  i,j  =  1,2,  and  Ai  is  real  analytic  in  (wi,w2)  except  at  (0,0); 
2 

(H2)  5Z  (II^IU-(O)  +  llu;j2^;'4«lli00(n))  <  +oo, 

i,j=l 

where  Q,  is  any  compact  set  containing  the  umbilic  point  P. 


Theorem  6.1.  Assume  that  (1.1)  satisfies  (HI)  and  (H2).  Suppose 
that  a  Young  measure  v  6  Proh.(R2)  satisfies 

(6.1)  supp  u(wx,w2)  C  [tuf,u;+]  x  [u;~, toj], 

(6.2)  <  v,  77^2  -  %5i  >=<  ><  v,q2  >  -  <  u,q2  ><  v,qx  >, 

for  all  C2  Goursat  entropy-entropy  pairs  (//,,  g,),  i  =  1, 2,  as  constructed 
in  Section  4  for  the  systems  (1.1)  with  an  isolated  umbilic  point  P. 
Then 

supp  v(wuw2)  C  {  (wuw2)  j  Xiw.  =  0,  i  =  1,2  }. 


Remark.  For  systems  with  more  than  one  isolated  umbilic  points,  a  variant  of 
Theorem  6.1  will  be  provided  in  [CK2]. 

Combining  Theorem  6.1  with  Theorem  5.3.1,  v.e  have 

Theorem  6.2.  Suppose  that  U((x,t )  are  measurable  functions  satis¬ 
fying 


(6.3)  ||*7%~  <  M  <  +oo, 

(6.4)  *l(U()t  +  q(U() x  compact  in  H^, 

for  any  C 2  Goursat  entropy-entropy  pair  ( r),q )  to  the  systems  (1.1) 
with  an  isolated  umbilic  point.  Suppose  that  (1.1)  satisfies  (Hi)  and 
(H2)  and  that  Atw.  ^  0,  i  —  1,2.  Then 


(6.5) 


(wi(Ue),w2(Uf))  — ►  (wi(x,t),w2(x,t)), 


a.e. 
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where  (w\{U),w2(U))  axe  the  Riemann  invariants. 


Remark.  Theorem  6.1  and  Theorem  6.2  provide  a  framework  for  the  compact¬ 
ness  of  the  sequence  of  the  Riemann  invariants,  determined  by  the  sequence  of 
approximate  solutions  satisfying  the  conditions  (6.3)  and  (6.4). 

We  sketch  the  ideas  in  the  proof  of  Theorem  6.1.  Without  loss  of  generality, 
we  assume  that  the  umbilic  point  P  =  (0,0).  Suppose  that  R  is  the  minimal 
rectangle  in  (uq,^)  space  containing  the  support  of  v.  There  are  two  cases: 

1.  R  does  not  contain  the  umbilic  point  (0,0); 

2.  R  contains  the  umbilic  point  (0, 0). 

We  will  treat  Case  2.  Case  1  is  similar  and  less  complicated.  We  also  assume 
that  R  is  not  a  line  segment  parallel  to  any  axis  (this  case  can  be  dealt  with  using 
the  method  shown  below),  thus,  R  =  [tu^O]  x  [0,  w^]  for  some  w <  0  <  u?^. 

Let  8  >  0  denotes  the  constant  used  in  the  construction  of  entropies  in 
Section  4.  Assume  that  it  is  chosen  a  priori  to  satisfy  w <  —  6  <  0. 

The  proof  will  rely  on  the  following  Propositions  and  Lemmas. 


Proposition  6.3.  Consider  the  entropy  function  constructed  in  Propo¬ 
sition  4.5.2.  Assume  that  w *  €  (tuf,—  8)  is  a  fixed  constant  and  that 
the  Goursat  data  6  is  supported  either  on  the  interval  ,w*)  (west 
type  entropy  with  limit  w*)  or  (it?*,  —8)  ( east  type  entropy  with  limit 
w*).  Then  the  entropy  rj  and  its  flux  q  admit  integral  representations 
of  the  form 


fWX 

r/(it?i,it?2)  =  l(w1,w2)0(wi)  +  /  J(wi,w2;y)  dy, 

Jw[ 

yu'i 

q(wi,w2)  =  K(wi,w2)8(wi)  +  /  I(u?i,u?2;y)  dy, 

J  wT 


where  /,  J ,  if,  L  are  smooth  as  long  as  y  <  0  and  w i  <  0.  Moreover, 

I>  1. 
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Proposition  6.4.  Let  (rpq)  be  of  type  east  with  limit  a  —  e1  and  let 
( fj,q )  be  of  type  west  with  limit  a  +  e' .  Suppose  that  for  all  e'  >  0 
sufficiently  small,  <  v,rj q  —  fjq  >=  0.  Then 

d\ 

supp  V  C  {  (u>i,u>2)  |  ~~(a,W2)  =  o  }. 


If  we  assume  also  that  u>2)  does  not  vanish.  Then  supp  u  D 

{(wi^W2)\w2  >  0,  w\  =  a}  is  empty. 

The  proof  of  Proposition  6.4  consists  of  a  construction  of  the  trace  of  u  on 
the  line  w\  =  a.  The  argument  is  quite  standard  and  we  refer  the  reader  to  [Se]. 

For  the  sake  of  simplicity,  we  assume  from  now  on  that  dwi  ^  i  —  1, 2,  is 
satisfied.  The  general  case  is  only  slightly  more  tedious.  With  this  assumption, 
v  can  be  reduced  to  a  point  mass  in  the  Riemann  plane  below. 

Next,  using  Proposition  6.3  and  some  properties  of  the  kernel  functions 
I,  J,  K ,  L ,  one  can  show  that 

Proposition  6.5.  Fix  6  >  0.  Let  w\  satisfy  w f  <  —  8  <  0.  If 

for  all  e^st  type  entropy  r\  with  limit  w*  we  have  <  v,r)  >=  0  then 
supp  v  fl  {(toi,ie2)|w2  >  0, re*  <  u>i  <  —6}  is  empty. 

Now,  we  will  make  use  of  Propositions  6.3  — 6.5  to  reduce  u  to  a  point  mass. 
The  reduction  process  will  take  several  steps.  First,  we  show  that  the  support  of 
v  must  concentrate  only  at  the  four  corners  of  R,  i.e.,  v  is  the  sum  of  four  delta 
functions.  Then  we  further  reduce  the  number  of  delta  functions  from  four  to 
one. 


It  turns  out  that  after  some  detailed  analysis,  Propositions  6.4  and  6.5  allow 
us  to  concentrate  the  support  of  v  on  the  extreme  left  vertical  edge  of  R  and  an 
arbitrarily  narrow  vertical  strip  containing  the  umbilic  point. 

Proposition  6.6.  supp  u  C  {(u;i,iy2)|0  <  iu2  <  w+,  and  uq  = 
wf  or  —  8  <  w\  <  0}. 
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By  using  entropies  of  type  north  and  south  and  by  performing  a  similar 
reduction  process  on  horizontal  line  segments  interior  to  the  rectangle  R ,  we 
conclude  that  suppu  is  concentrated  on  a  horizontal  strip  of  width  6  and  cn 
the  line  uq  =  w* ■  Due  to  the  complete  similarity  of  the  proof,  we  omit  this 
argument  but  summarize  the  result  as  follows: 

Proposition  6.7.  supp  v  C  {(uq,u;2)|0  <  w\  <  ,  and  u>2  = 

w£  or  8 >  W2  >0}. 

Proposition  6.6  and  6.7  imply  that  the  support  of  v  is  in  fact  concentrated  on 
the  square  {(uq, u>2)|0  <  uq  <  8,  —  8  <  uq  <  0},  and  the  points  (w^~, 0), (0,  u^), 
and  Now,  8  >  0  is  small  but  arbitrarily  fixed.  Let  8  — *  0.  We  obtain 

Corollary  6.8.  The  support  of  v  is  concentrated  at  the  four  comers 

of  R,  (0, 0),  (0,  ),  0),  and  (w^ )■ 

Knowing  that  the  Young  measure  v  is  a  sum  of  four  delta  functions  at  the 
corners  of  R ,  we  reduce  v  further.  Let  A,,  i  =  1,2, 3, 4,  denote  the  corners  of 
R,  Ai  =  (0,u^),A2  =  (0,0), A3  =  (u>f,0),A4  = 

By  Corollary  6.8, 

^  M Ah 

l<i<4 

where 

Pi  >  0,  i  =  1,2, 3,4,  and  ^  /?,  =  1. 

1<»<4 

To  reduce  v  further,  we  apply  an  variant  of  Theorem  6.1  in  [Se]. 

Proposition  6.9.  Suppose  that  Pi  >  0  for  all  i  =  1,2, 3,4.  Then 
dWlX1  =  0,  on  the  line  segments  A\A 4,  A2A3; 
dw2X.2  =  0,  on  the  line  segments  A1A2,  A3A4. 

Now,  by  assumption,  dWl\\  and  dW2\2  are  nonzero  everywhere. 

This  contradicts  the  conclusion  of  the  above  Proposition.  Thus  we 

have 
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Lemma  6.10.  There  exists  an  i,  1  <  i  <  4,  such  that  fa  =  0. 


Therefore,  v  is  only  a  sum  of  at  most  three  Dirac  masses.  Further  reduction 
of  v  can  be  achieved  by  combining  polynomial  entropies  (Proposition  3.1)  and 
Goursat  entropies.  Finally,  we  conclude  that  the  Young  measure  u  is  a  point 
mass  on  the  (nq,  u>2)-plane. 

This  completes  the  sketch  of  the  proof  of  Theorem  6.1. 


7.  Applications  of  the  Compactness  Framework 

In  this  section,  we  apply  the  compactness  framework  theorem  to  the  qua¬ 
dratic  flux  system  in  region  IV  (A  >  0).  We  prove  the  strong  convergence  of 
approximate  solution  sequences  constructed  by  the  viscosity  method,  the  Go¬ 
dunov  scheme,  and  the  Lax  Friedrichs  scheme  to  the  weak  entropy  solutions  of 
(2.5)  with  large  initial  data  for  all  positive  time.  As  a  corollary,  we  obtain  the 
global  existence  of  weak  solutions  to  the  Cauchy  problem  of  (2.5)  with  large 
data  in  L°°. 

We  first  verify  L°°  apriori  estimates  for  these  approximate  solution  se¬ 
quences.  We  then  apply  the  compactness  framework  (Theorem  6.2)  together 
with  Theorems  4.5.2,  4.5.3,  and  4.5.6  on  the  existence  of  regular  entropies.  To 
this  end,  the  method  of  invariant  regions  will  be  used  for  the  viscosity  approxi¬ 
mate  sequence.  In  the  case  of  the  finite  difference  approximations,  it  is  necessary 
to  establish  such  apriori  estimates  for  general  Riemann  problem  solutions  to  (2.5) 
which  are  the  building  blocks  of  these  difference  schemes.  We  remark  that  the 
Riemann  problem  for  (2.5)  has  not  been  solved  by  purely  analytical  methods 
before.  Here  we  use  the  viscosity  method  to  obtain  existence  of  solutions  to  the 
Riemann  problem  on  arbitrary  compact  sets  in  the  (x,<)-plane.  At  the  same 
time,  we  obtain  the  desired  L°°  apriori  estimates. 


7.1  Convergence  of  the  Viscosity  Method 
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Consider  the  Cauchy  problem 

dtu  +  dx(dC(U))  =  o, 

U(x,0)  =  Uo(x), 

C(U)  =  ^(^au3  +  bu2v  +  uv2), 

L  O 

(7.1)  A  >  0. 


Using  the  convexity  properties  of  the  Rj  curves  or  the  explicit  form  of  Wj, 
we  obtain 

Proposition  7.1.1.  Consider  the  quadratic  flux  system  (2.5).  Assume 
that  we  are  in  region  TV,  i.e.,  A  >  0.  Then  for  any  pair  of  constants, 

Ci  <  0  <  Ci,  QCic2  and  QClCs  Hit,  with  Ik  as  defined  in  Proposi¬ 
tion  5.1.1,  axe  convex  invariant  regions  for  the  parabolic  systems  with 
viscosity  associated  with  (2.5). 

Using  Proposition  7.1.1,  we  obtain  an  L°°  apriori  bound  for  viscosity  ap¬ 
proximate  solutions  Uf . 

Proposition  7.1.2.  Consider  the  Cauchy  problem  (7.1)  with  Cauchy 
data  Uo  €  Tk,  and  its  associated  viscosity  approximation  {Uf}t>o. 
Suppose  that  there  is  a  constant  U  such  that  Uq—U  6  fi2(R)nL°°(R). 
Then,  for  any  e  >  0,  Uf(x,  t )  is  well-defined  for  all  (x,  t)  and  moreover, 
|^(x,t)|<||U0|Uco. 

Combining  Proposition  7.1.2,  Theorem  6.2,  and  Theorem  4.5.6, 
we  obtain 

Theorem  7.1.3.  Consider  the  Cauchy  problem  (7.1)  with  Cauchy 
data  Uo  €  Ijt,  and  its  associated  viscosity  approximation  {U(  }e>o- 
Suppose  that  A  >  0.  Suppose  that  there  is  a  constant  U  such  that 
Uo  —  U  €  L2{ R)  fl  ir°°(R).  Then,  as  e  —>  0+,  there  exists  a  subse¬ 
quence  ofUt(x,t)  that  converges  a.e.  (x,  t)  to  a  global  weak  solution 
of  the  hyperbolic  system  (7.1). 
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7.2  Convergence  of  Finite  Difference  Schemes 

The  Godunov  and  the  Lax- Friedrichs  schemes  use  the  solutions  of  the  Rie- 
mann  problem  as  their  building  blocks.  In  the  case  of  the  quadratic  flux  system 
(2.5),  explicit  constructions  of  the  Riemann  problem  solution  prove  to  be  ex¬ 
tremely  complicated.  And  the  Riemann  problem  has  never  been  solved  before 
by  pure  analytical  methods.  In  order  to  prove  the  convergence  of  these  finite  dif¬ 
ference  schemes  for  the  quadratic  flux  system  with  arbitrary  data  in  L°°,  we  first 
establish  the  existence  of  solutions  for  Riemann  solutions  and  their  L°°  apriori 
bounds.  To  this  end,  we  apply  the  viscosity  method  studied  in  Subsection  7.1. 

Proposition  7.2.1.  Consider  the  Riemann  problem  for  (2.5): 
dtU  +  dx(dC(U ))  =  0, 

(7.2)  C(U)  =  ^(^au3  +  bu2v  +  uv2),  a/l-ffe2, 

z  o 

with 


(7.3) 


Ul,  x  <  0, 
Ur,  x  >  0. 


Here,  Ul,Ur  €  Ik  are  constants.  Suppose  that  A  >  0.  Let  II  C 
{  (x,t)  |  x  €  R,  t  >  0  }  be  a  compact  set.  Then  there  exists  a  weak 
solution  to  the  Riemann  problem  (7.2)  on  the  domain  II. 


The  proof  of  Proposition  7.2.1  makes  use  of  the  finiteness  of  propagation 
speed  for  solutions  of  (7.2)  and  (7.3)  and  of  Theorem  7.1.3.  The  invariant 
regions  for  the  viscosity  approximate  solutions  are  convex.  Thus,  the  L°°  bound 
that  the  Riemann  solutions  inherit  from  the  viscosity  method  axe  preserved 
during  the  averaging  processes  in  the  intermediate  steps  in  the  finite  difference 
approximations. 


Proposition  7.2.2.  Consider  the  Cauchy  problem  (7.1)  with  Cauchy 
data  Uo  6  Ik  for  the  quadratic  flux  system.  Assume  that  Uo  G  L°°. 
Let  {[/*}*>o  be  a  sequence  of  approximate  solutions  to  (7.1)  generated 
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by  the  Lax-Friedrichs  scheme  or  the  Godunov  scheme.  Then 
is  well-defined  and  \Ui(x,t)\  <  ||^o|U°°- 

We  obtain,  after  combining  the  above  propositions  with  Theorem  4.5.6  and 
the  analysis  in  Section  6, 

Theorem  7.2.3.  Consider  the  Cauchy  problem  (7.1)  for  the  quadratic 
fiux  system  with  Cauchy  data  Uq  6  2*.  Assume  that  Uq  €  L°°.  Let 
{Ue}t>0  be  a  sequence  of  approximate  solutions  to  (7.1)  generated  by 
the  Lax-Friedrichs  scheme  or  the  Godunov  scheme.  Then,  as  i  — ►  0+, 
there  exists  a  subsequence  of  that  converges  to  a  global  weak 

solution  of  (7.1)  a.e.  in  ( x,t ). 

Remark.  Our  method  yields  the  same  convergence  and  existence  theorems  for 
another  nonstrictly  hyperbolic  system  (a  special  case  of  such  system  was  studied 
in  [Rb]):  (1.1)  with  F  =  (au2  uv)T,  a  >  1,  for  all  even  regular,  superlinear 

function  ^p(v)  satisfying  ip'(v)v  >  0  and  (p"(v)v  >  0  when  v  ^  0. 
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